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Abstract. Empirical star formation laws from the last 20 years are re- 
viewed with a comparison to simulations. The current form in main 
galaxy disks has a linear relationship between the star formation rate 
per unit area and the molecular cloud mass per unit area with a 
timescale for molecular gas conversion of about 2 Gyr. The local ratio of 
molecular mass to atomic mass scales nearly linearly with pressure, as 
determined from the weight of the gas layer in the galaxy. In the outer 
parts of galaxies and in dwarf irregular galaxies, the disk can be dom- 
inated by atomic hydrogen and the star formation rate per unit area 
becomes directly proportional to the total gas mass per unit area, with 
a consumption time of about 100 Gyr. The importance of a threshold 
for gravitational instabilities is not clear. Observations suggest such 
a threshold is not always important, while simulations generally show 
that it is. The threshold is difficult to evaluate because it is sensitive 
to magnetic and viscous forces, the presence of spiral waves and other 
local effects, and the equation of state. 



1 Introduction: the Kennicutt-Schmidt law of Star Formation 



• • , Star formation and stellar evolution are such important drivers of galactic evolution 
m ^ ' that empirical laws to determine the star formation rate have been investigated 
, for over 50 years. The results have never been very precise because star formation 
!_h ' spans a wide range of scales, from cluster-forming cores to molecular clouds to the 
5^ i whole interstellar medium. 

On the scale of a galaxy, the first idea was a proposed con nection between 
the total star formation rate and the mass of interstellar gas. ISchmidtl (Il959h 



derived the star formation rate (SFR) over the history of the Milky Way assuming 
a constant initial luminosity function for stars, $(My), a stellar lifetime function 
T(My), a gas return per star equal to all of the stellar mass above 0.7 M , and a 
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star formation rate f(t) that scales with a power n of the gas mass, Ma{t). Then 
/(t)S MV , I , (Mv) = C[M G (t)] n , for a summation E MV over all stellar types. 

Schmidt gave analytical solutions for n = 0, 1, 2. He noted that a scale height 
for HI of 144 pc, a scale height for Cepheids of 80 pc, and a scale height for clusters 
of 58 pc gave n = 2 to 3. The white dwarf count gave n > 2, the He abundance 
suggested n — 2, the uniformity of HI suggested n > 2, and the cluster mass 
function gave n — 1 to 2. Schmidt also suggested that with n = 2, dense galaxies 
like ellipticals should now have less gas than low-density galaxies like the LMC. His 
final comment was "It is hoped to study the evolution of galaxies in more detail 
in the future." Following ISchmidtl ( 1959T) . many authors derived scaling relations 
between the average s urface density of star formation, Esfr, and the average 
surface density of gas. Buat. Deharveng fc Donas! ( 19891 ) included molecular and 
atomic gas and determined star formation rates from the UV flux corrected for 
Milky Way an d internal extinction. They assumed a constant H2/CO ratio and 



Scald (jl986l ) IMF. The result was a good correlation between the average star 



formation rate in a sample of 28 galaxies and the 1.65 ± .16 po wer of the average 
total gas surface density. In the same year, Kennicuttl (1989) used Ha for star 
formation, and HI and CO for the gas with a constant H2/CO conversion factor, 
and determined star formation rates both as a function of galactocentric radius 
and averaged over whole galaxy disks. For whole galaxies, the average Ha flux 
scaled with the average gas surface density to a power between 1 and 2; there 
was a lot of scatter in this relation and the correlation was better for HI than H2 . 
More interesting was Kennicutt's (1 989) resu l t tha t the star formation rate had 
an abrupt cutoff in radius where the iToomrd (|l964l ) stability condition indicated 
the onset of gravitationally stable gas. Kennicutt derived a threshold gas column 
density for star formation, E cr i t = acrre/(3.36G) for a = 0.7; a is the velocity 
dispersion of the gas; k is the epicyclic frequency, and G is the gravitational 

constant. 

In a second study. iKennicutt ( 19981 ) examined the disk-average star formation 
rates using a larger sample of galaxies with Ha, HI, and CO. He found that for 
normal galaxies, the slope of the SFR-surface density relation ranged between 
1.3 to 2.5, depending on how the slope was measured; there was a lot of scatter. 
When starburst galaxies with molecular surface densities in excess of 100 M & were 
included, the overall slope became better defined and was around 1.4. This paper 
also found a good correlation with a star formation rate that scaled directly with 
the average surface density of gas and inversely with the rotation period of the 
dis k. This second law s uggested that large-scale dynamical processes are involved. 
Hun ter et al. (|l998l) considered the same type of analysis fo r dwarf Irr e gulars 



and derived a critical surface density that was lower than the Kennicuttl ( 19891) 
value by a factor of ~ 2. This meant that stars form in more stable gas in dwarf 
irr egulars compar e d to s pirals . 

iBoissiet et al. (|2003l) compared Ssfr and S gas versus radius in 16 resolved 
galaxies with three theoretical expressions. The best fits were a SFR depen- 
dence on the gas su r face d ensity as Xsfr oc Sg a g 6 , a more dynamical law from 
iBoissier fc Prantzosl (|l999h which gave the fit Ssfr °c f° r rotation 
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speed V and radius R, and a third type o f law fr om Dopita fc Rvderl (1994), which 
fit to E SFR oc E^f/Eg-f 1 . iBoissiet et al.l (120031) assumed that H 2 /CO varied with 
radius as the metallicity (jBoselli et all 120021) . Their conclusion was that the three 
laws a re equally good, and that for the pure gas law, n > 1.4. IBoissiet et al 



(2003) also looked for a star formation threshold in the Milky Way. They deter 



mined E/E cr i t using both pure-gas for E cr it and a gas+star E cr it from I Wang fe Silk 



( 19941) . They found that the gas+star E cr i t gave the best threshold for determin- 



ing where star formation occurs. The gas alone was sub-threshold throughout the 

disk. 

Zasov fc Smirnova (2005) showed that a threshold like E crit may be used to 



determine the gas fraction in galaxies. If all galaxies have E(HI) approximately at 
the critical E cr i t = ana/irG, which is proportional to V/R from k, then M gas = 
J R 2TTRE a ndR oc VR. This was shown to be the case from observations. They 
also considered that the total mass is M tot oc V 2 R, in which case M tot /M gas oc V, 
the rotation speed. This was also shown to be confirmed by observations. In 
their interpretation, small galaxies are more gas-rich than large galaxies because 
all galaxies have their g as column densit i es clo se to the surface density threshold. 
For the Milky Way. iMisiriotis et al.l (|2006l ) used COBE/DIRBE observations 



to get both the gas and dust distributions and the SFR distribution. They found a 
gas-law slope of 2.18 ±0.20, which they claimed w as similar to Kenni cutt's (1998) 
bivariate fit slope n = 2.5 for normal galaxies. iLuna et al. (200G) determined 



the Milky Way SFR from IRAS point sources and the CO surface density from a 
southern hemisphere survey (assuming constant H2/CO). They found star forma- 
tion concentrated in low-shear spiral arms and sug gested an additio nal dependence 
on shear. Overall they derived Esfr ~ E ga 2±0 - 2 . IVorobvovl ( 2003h also suggested 



a shear dependence for the SFR based on observations of the Cartwheel galaxy, 
where there is an inner ring of star formation with high shear that is too faint for 
the normal Kennicutt law, given the gas column density. 

1.1 TheQ Threshold 

A threshold for gravitational instabilities in rotating disks has been derived for 
various ideal cases. For an infinitely thin disk of isothermal gas, the dispersion 
relation for radial waves is cj 2 = fc 2 er 2 — 2-7rGEfc + K 2 . Solving for the fastest growth 
rate uj gives the wavenumber at peak growth, k = irGY^/a 2 , and the wavelength, 
A = 2er 2 /GE, which is on the order of a kiloparsec in main galaxy disks. The 
dominant unstable mass is M ~ (A/2) 2 E = er 4 /G 2 E ~ 10 7 Mq in local spirals. 
The peak rate is given by 

"peak = -(ttGE/ cr 2 ) 2 + k 2 = -(^GE/a 2 ) 2 (l - Q 2 ) (1.1) 

which requires Q = na/irGY, < 1 for instability (i.e., when w 2 eafe < 0). 

Disk thickness weakens the gravitational force in the in-plane direction by an 
amount that depen ds on wavenumber, approximately as l/(l + fcg) for exponential 
scale height H (e.g.. lElmegreenl . ll987l : lkim fc Ostrikerll2007l ). Typically, k~l/H, 



so this weakening can slow the instability by a factor of ~ 2, and it can make 
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the disk slightly more stable by a factor of 2 in Q. On the other hand, cooling 
during condensation decreases the effective value of the velocity dispersion, which 
should really be written for adiabatic index 7 that appears in the relation 

SP oc Sp 1 with pressure P and density p. If P is nearly constant for changes in p, 
as often observed, then 7 ~ 0. Mversl (1978) found 7 ~ 0.25 for various thermal 
temperatures at interstellar densities between 0.1 cm -3 and 100 cm -3 . Thus the 
effects of disk thickness and a soft equation of state partially compensate for each 
other. 

There is also a Q threshold for the collapse of an expanding shell of gas 
( Elmegreen. Palous fc Ehleroval 120021 ). Pressures from OB associations form giant 
shells of gas and cause them to expand. Eventually they go unstable when the 
accumulated gas is cold and massive enough, provided the induced rotation and 
shear from Coriolis forces are small. Considering thousands of initial conditions, 
these authors found that a sensitive indicator of whether collapse occurs before 
the shell disperses is the value of Q in the local galaxy disk, i.e., independent of 
the shell itself. The fraction / of shells that collapsed scaled inversely with Q as 
/~0.5-0.41og 10 Q. 

The Toomre Q parameter is also likely to play a role in the occurrence o f 
instabilities in turbulence-compressed gas on a galactic scale ( ElmegreenL 2002 ). 
Isothermal compression has to include a mass comparable to the ambient Jeans 
mass, Mjoans, in order to trigger instabilities. The turbulent outer scale in the 
galaxy is comparable to the Jeans length, £j G ans, which is about the galactic gas 
scale height, H . If the compression distance exceeds the epicyclic length, then 
Coriolis forces spin up the compressed gas, leading to resistance from centrifugal 
forces. So instability needs Lj ean s < -^epicycle, which means Q < 1, since £j ca ns ~ 
H — a 2 /7rGS. The epicyclic length is L op icycie ~ cr/k, so £j Gans /L p icyc i = Q. 

The dimensionless parameter Q measures the ratio of the centrifugal force from 
the Coriolis spin-up of a condensing gas perturbation to the self-gravitational force, 
on the scale where gravity and pressure forces are equal, which is the Jeans length. 
The derivation of Q assumes that angular momentum is conserved, so the Cori- 
olis force spins up the gas to the maximum possible extent. When Q > 1, a 
condensing perturbation on the scale of the Jeans length spins up so fast that 
its centrifugal force pulls it apart against self-gravity. Larger-scale perturbations 
have the same self-gravitational acceleration (which scales with E) and stronger 
Coriolis acceleration (which scales with K 2 /k); smaller-scale perturbations have 
stronger accelerations from pressure. If angular momentum is not conserved, then 
the disk can be unstable for a wider range of Q because there is less spin up 
during condensation. For example, the Coriolis force can be resisted by magnetic 
tension or viscosity and then the angular momentum in a cond ensing cloud will 
get stripped away. Th i s removes the Q threshold completely dChandrasekhar , 
1954 IStephensonl . Il96lt iLvnden-Bell Il966t iHunter fc Horakl |l983|). In the mag- 



netic case, the result is the Magneto- Jeans instability, which can dominate the 
gas condensation i n low-shear environments like spiral arms and some inner disks 
(|Elmegreenl . Il987l Il99ll Il994t iKim fc OstrikerL I2001 120021 iKim et all 12002ft . For 
the viscous case. iGammid (|l996| ) showed that for Q close to but larger than 1, i.e., 
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in the otherwise stable regime, viscosity can make the gas unstable with a growth 
rate equal to nearly one-third of the full rate for a normally unstable (Q < 1) disk. 
A dimensionless p arameter for viscosity v is vk 3 jG 2 Y? , which is ~ 11 according 
to iGammid (|l996l) . This is a large value indicating that galaxy gas disks should 
be destabilized by viscosity. An important dimensionless parameter for magnetic 
tension is B 2 / (ttGT, 2 ) ~ 8, which is also large enough to be important. Thus gas 
disks should be generally unstable to form small spiral arms and clouds, even with 
moderately stable Q, although the growth rate can be low if Q is large. 



1.2 Modern Versions of the KS Law with ~ 1.5 slope 



Kennicutt et al.l ([2007Q studied the local star formation law in M51 with 0.5-2 kpc 
resolution using Pa-a and 24/i+Ha lines for the SFR, and a constant conversion 
factor for CO to H 2 . There was a correlation, mostly from the radial variation 
of both SFR and gas surface density, with a slope of 1.56 ± 0.04. There was 
no correlation with S(HI) alone, as this atomic component had about constant 
column density (~ 10 Mq). The correlation with molecules alone was about the 
same as the total gas correlation. 



Lerov et al. (120051) studied dwarf galaxies and found that they have a molec- 



ular KS index of 1.3 ± 0.1, indistinguishable from that of spirals, except with a 
continuation to lower central H 2 column densities (i.e., down to ~ 10 Mq pc~ 2 ). 



Hever et al.l (|2004l ) found a slope n — 1.36 for Hsfr versus £(H 2 ) in M33, 
where the molecular fraction, ,f m oi is small. The co rrelation with the total gas 
was much steeper. More recently, Verlev et al. ( 2010f ) studied M33 again and got 



Ssfr oc £g 2 for n = 1 to 2, and Esfr oc £™ otal gas for n = 2 to 4. The steepening 
for total gas is again because Shi is about constant, so the slope from HI alone 
is nearly infinite. This correlation is dominated by the radial variations in both 
quantities, as it is a point-by-point evaluation throughout the disk. Radial changes 
in metallicit y, spiral arm activa tion, tidal density, and so on, are part of the total 
correlation. Verlev et al.l (|2010h also try other laws, such as Ssfr oc (Sh2Pism)™, 



for which n = 1.16±0.04, and £ S fr oc pf SM , for which n = 1.07±0.02. These differ 
by considering the conversion from column density to midplane density, using a 
derivation of the gaseous scale height. The first of these would have a slope of unity 
if the star formation rate per unit molecular gas mass were proportional to the 
dynamical rate at the average local (total) gas density. The second has the form 
of the original Schmidt law, which depends only on density. To remove possible 
effects of CO to H2 conversion, Verley et al. also looked for a spatial correlation 
with the 160 /x opacity, rigo, which is a measure of the total gas column density 
independent of molecule formation. They found Ssfr oc t™ 60 for n = 1.13 ± 0.02, 
although the correlation was not a single power law but a 2-component power law 
with a shallow part (slope ~ 0.5) at low opacity (ri6o < 10 -4 ) and a steep part 
(slope ~ 2) at high opacity. 
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1.3 Explanations for the 1.5 slope 



Prior to around 2008, the popular form of the KS law had a slope of around 1.5 
when Ssfr was plotted versus total gas column density on a log-log scale. This 
follows from a dynamical model of star formation in which the SFR per unit area 
equals the available gas mass per unit area multiplied by the rate at which this 
gas mass gets converted into stars, taken to be the dynamical rate, 



SsFR ~ eEgas (GjOgas) 



1/2 



(1.2) 



If the gas scale height is constant, then S gas oc p gas and Ssfr oc S ga g. In the model 
of star formation where star-forming clouds are made by large-scale gravitational 
instabilities, this 1.5 power law would work only where the Toomre instability 
condition, Q < 1.4, is satisfied. Such a model accounts for the Kennicutt (1989, 
1998) law with the Q < 1.4 threshold. 

Several computer simulations have shown this dynamical effect. Li et al.l ( 20061) 
did SPH simulations of galaxy disks with self-gravity forming sink particles at den- 
sities larger than 10 3 cm -3 . They found a Q threshold for sink particle formation , 
and had a nice fit to the KS law with a slope of ~ 1.5. iTasker fc Brvanl (|2006l) 
ran ENZO, a 3D adaptive mesh code, with star formation at various efficiencies, 
various temperature floors in the cooling function, and various threshold densities. 
Some models had a low efficiency with a low threshold density and other models 
had a high efficiency with a high threshold density. Some of their models had 
feedback from young stars. They also got a KS slope of ~ 1. 5 for both global and 
local star formation, regardless of the details in the models. iKravtsov (|2003f) did 
cosmological simulations using N-body techniques in an Eulerian adaptive mesh. 
He assumed a constant efficiency of star formation at high gas density, and star 
formation only in the densest regions (n > 50 cm -3 , the resolution li mit), which 
are in the tail of the density prob a bility distr i bution function (pdf; cf. lElmegreen . 
20021 : iKrumholz fc McKeeL 120051 ). iKravtsovl (120031) go t the KS law with a slope 
of 1.4 for total gas surface density. Wada fc Norman ( 20071) did a similar thing, 
using the fraction of the mass at a density greater than a critical value from the 
pdf (/Ocrit = 10 3 cm -3 ) t o determine the star formatio n rate. Their analytical 



result had a slope of 1.5. Harfst. Theis fc Hensler ( 20061 ) had a code with a hier- 
archical tree for tracking interacting star particles, SPH for the diffuse gas, and 
sticky particles for the clouds. They included mass exchange by condensation and 
evaporation, mass exchange from stars to clouds (via PNe) and from stars to dif- 
fuse gas (SNe), and from clouds into stars during star formation. New clouds were 
formed in expanding shells. Their KS slope was 1.7 ±0.1. They also got a drop in 
Ssfr at low S gas , not fr om a Q threshold but from an inability of the gas to cool 
and form a thin disk (cf. iBurkert et al. . 1992 : Elmegreen fc Parravano , 1994 ). 



2 The Molecular Star Formation Law 



The star formatio n law may also be writte n as a linear rela t ion fo r molecules, with 



Ssfr oc S^j 2 fe.g.. iRownd fc Youngl . fl999l ) . IWong fc Blital (|2002t ) found a SFR in 
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direct proportion to molecular cloud density (n — 1), and suggested that the n = 
1.4 KS law came from changes in the molecular fraction, f mo \ — £112/ (^hi + ^m)- 
They assumed that H 2 /CO was constant and determined the combined index n' = 
n mo\ (1 + din /mol/^hiSgas) where n mo \ = 1 and / mo i increases with pressure, P. 
They measured d1nf mo \/d1nP ~ 0.2, and if P oc £ 2 as , then ciln/ mo i/cilnE gas ~ 
0.4. This gives the KS n = 1.4 law for total gas. Wong & Blitz also suggested 
that the stability parameter Q was not a good threshold for star formation, but a 
better mea sure of the gas fracti o n in t he sense that a high Q corresponds to a low 
/Etot- iBlitz fc Rosolowskil (|2006l ) showed for a wider sample of 13 galaxies 
that the molecular ratio, i? mo i = Ejk/Shi, sca l es about linearly with the total 
ISM pressure. Interacting galaxies had slightly higher i? mo i for a given P, but 
among interacting galaxies, the correlation was still present. 

A large study of HI, CO, and star f ormation rat e s from GA LEX ultraviolet an d 
Spitzer 24// observations was made bv lBigiel et al. ( 20081 ) and Lerov et al. ( 2008 ). 
They considered the local star formation law with a resolution of 750 pc. Bigiel 
et al. found that Esfr oc £co> an d tha t the timescale for conversion from H 2 to 
stars was about 2 Gyr. Figure Q] (from Bigiel et al. . 20081) shows an example of 
how much better the SFR scales with CO than either HI or the total gas. The CO 
and SFR maps of NGC 6946 resemble each other closely, and neither resembles 
the HI map. Bigiel et al. also found that Ehi saturates to ~ 9 M Q pc~ 2 . When 
plotting SgFR over a wide range of £hi+H2, they found a slope of unity in the 
molecular range, Shi+h2 > 9 M Q pc -2 , and higher slope in the atomic range 
(Shi+H2 < 9 Mq pc -2 ). Figure [2] shows the summed distribution of SFR per unit 
area versus total gas column density in 7 spiral galaxies. There is a linear part at 
high column density and a steeper part at low column density. 

Dwarf galaxies look like the outer parts of spirals in the Bigiel et al. survey, 
occupying the steeper part of the SgpR — S gas diagram at low E gas . At higher 
Shi+H2, the survey did not have new data, but Bigiel et al. suggested, based 




Fig. 1. Maps of HI, CO and SFR in NGC 6946 with HI on t he left, CO in the m iddle 
and SFR on the right, all convolved to 750 pc resolution (from Bigiel et all (|2008l ). The 
circle is the optical radius at 25 mag arcsec -2 . 
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on Kennicutt's (1998) starburst result, that perhaps the KS law turned up to a 
steeper slope (n ~ 1.4) in a third regime of star formation where £h2 exceeds the 
standard column density of a single molecular cloud (around 100 M pc~ 2 ). 



Leroy et al.l ( 20081 ) compared these new survey results to various theoretical 



models. They found that the star formation time in CO-rich gas is universally 
1.9 Gyr, independent of the average local free fall or orbital time, the midplane 
gas pressure, the state of gravitational stability of the disk with or without the 
inclusion of stars in the stability condition, and regardless of the rate of shear 
or the ability of a cold gas phase to form. Star formation depends only on the 
presence of molecules and it proceeds at a fixed rate per molecule. Leroy et al. also 
found that dwarf galaxies are forming stars at their average historical rate, whereas 
spirals are forming stars at about half of their average rate. In the outer disk, the 
SFR in HI drops with radius faster than the free fall time, suggesting self-gravity 
is not the lone driver. Also important are the phase balance between HI and H2, 
giant molecular cloud (GMC) destruction, stellar feedback, and other processes. 
These processes govern the presence of GMCs with an apparently constant star 
formation efficiency in each GMC. 

Unlike the star formation rate per molecule, the molecule-to-atom ratio does 
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Fig. 2. The distribution of SFR per unit area versus total gas column density, convolved 
to 750 pc, for 7 spiral galaxies (from[Bigicl et all l|200St ). There is a change in the slope 
from ~ 4 at £m+H2 < 9 Mq pc -2 (the vertical dashed line) in the outer disk to ~ 1 
at higher £m+H2 in the inner disk. The short-dashed lines correspond to gas depletion 
times of 0.1 Gyr, 1 Gyr, and 10 Gyr, from top to bottom. 
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correlate well with environmental parameters. iLerov et al. I (120081) showed approx- 
imately linear correlations with stellar surface density and interstellar pressure, an 
inverse squared dependence on the orbit time, and an exponential dependence on 
the galactic radius, like the rest of the disk, and with a comparable radial scale 
length. The molecular fraction is a smooth function of environmental parame- 
ters (e.g., pressure); no thresholds were seen. Disks seem to be marginally stable 
throughout. 

Leroy et al. concluded by noting that the HI-H2 transition in spirals typ- 
ically occurs at 0.43 ± 0.18 R25, which is about the same as where S sta rs = 
81 ± 25 Mq pc~ 2 , S gas = 14 ± 6 M pc~ 2 , P = 2.3 ± 1.5 x 10 4 /c B K cm" 3 , 
and Tot-bit = 1-8 ± 0.4 Gyr. There should be unobserved H 2 in dwarfs, according 
to the high star formation rates and low CO emissions there; in fact Leroy et al. 
estimate for dwarfs Ejj2 ~ 2Shi m the inner regions. 

Where Shi > Sh2, the star formation efficiency is proportional to S sta rs, mak- 
ing 

S SF r ~ S gas x [S stars /81 Mq pc- 2 ] /1.9 Gyr. (2.1) 



2.1 Theoretical Models for the Bigiel-Leroy Observations 



Krumholz et al.l ((2008) considered the molecule formation problem by starting 
with the radiative transfer of H2-dissociating radiation: dF/dz = —n<JdF—f}nn 2 R/fdh 
Here, F is the flux in Lyman- Werner bands that dissociate H2, n is the density 
(~ 30 cm -3 near the H2 transition), Cd is the dust cross section per H (10 -21 
cm 2 ), /hi is the fraction of n that is HI, R is the rate coefficient for formation of 



ISS 



H2 on grains (~ 3 x 10 cm 3 s , i.e., the formation rate is fmn R), and 
is the fraction of uv photon absorptions that dissociate H2 (~ 0.1). 

The solution to this radiative transfer equation is F(t) = (e~[ r ~ THI l — l) j\ 
where thi = ln(l + x); X = j dissO 'dcEg / 'nR is the ratio of absorption in dust to H2, 
and Eq is the free space photon number density (~ 7.5 x 10~ 4 cm -3 ). Krumholz et 
al. (2008) assume that the cold neutral density comes from two phase equilibrium. 
Then n scales with Eq and \ becomes nearly constant. From this they get the 
extinction, Ay, to the HI/H2 transition, the HI column density, Shi, and the 
molecular fraction in spherical cloud complexes as a function of the complex total 
column density. They do this also as a function of metallicity. 

After considering a galactic cloud population. iKrumholz et al~ ( 2009a ) derive 



Shi versus S to tai gas for different metallicities, and compare this with observations 
of galaxies having those metallicities. They do the same for H2. They also com- 
pare the observed versus the predicted correlation between H 2 /HI and pressure 
P. To do this, they use the observed S to tai gas and metallicity, and then compute 
i?H2 = Sho/Sht from the ory. This is plotted ver sus the observed pressure from 



Blitz fc Rosolowskl (|2006t ) and Leroy et al. ( 20081 ). The agreement is good 



Krumholz et al. (20091)) considered the star formation law, 



SsFR = ^IgasfmSFRff/tff 



(2.2) 
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where the star formation rate in a free fall time is the fraction of the gas that 
turns into stars in a free fall time, SFRg, divided by the free fall time, tg. This is 
SFR s /t s = (M 6 " a33 /0.8 Gyr) x Max[l, S gas /85 M pc" 2 ]; M 6 is the cloud mass 
in units of 10 6 M Q . This equation assumes that stars form in the high density 
tail of a log-normal density pdf, with the tail width given by the Mach number; a 
fraction of 0.3 of the dense gas mass goes into stars. The clouds are virialized and 
at uniform pressure until the galactic E gas exceeds the column density of a single 
GMC; then the pressure equals the galactic pressure. Also, the cloud complex 
mass is taken to be M$ = 37S gas /(85 Mq pc~ 2 ) from the galaxy Jeans mass. 

This theory for molecule form ation and sta r forma tion i n a galactic e n viron ment 
fits well to the observations by iBigiel et al.l ( 2008 ) and iLerov et al. ( 2008 ). It 
reproduces the low column density regime by having a low ratio of molecules 
to atoms at low pressure, it reproduces the intermediate column density regime 
by having a hxed star formation rate per molecule and an areal average star 
formation rate from the areal density of molecular clouds at constant pressure, 
and it reproduces the high column density regime by increasing the interstellar 
pressure, which makes the cloud density go up and the free fall time go down. A 
key point in their model is that molecular cloud pressures are constant in normal 
galaxy disks because they are set by HII region pressures (feedback) and not the 
galactic environment. In this sense, all GMCs have to be parts of shells or other 
active disturbances formed by high pressures. 

We k now that mole cular cloud pressures in the Milky Way are about constant 
from the iLarsonl (|198lh laws, which require this for virial ized clouds, but we don't 



really know the reason for it. It could be feedback, as iKrumholz et al. (2009b) 



suggest, or it could be the weight of t he HI shielding layer, which has a regulatory 
effect on pressure (jElmegreenl . 1989). This regulatory effect works because at 
high ambient pressure, the atomic density on the periphery of molecular clouds is 
high and so the required surface column density for H2 line self-shielding is low, 
and vice versa. The pressure at the bottom of the shielding layer, which is the 
molecular cloud surface pressure, scales directly with the column density of the 
shielding layer. Thus a lower intercloud pressure is compensated by a higher HI 
column density at the molecular cloud surface, making the molecular cloud surface 
pr essure somewhat uniform . 



Robertson &: Kravtsovl (|2008l ) simulated star formation in galaxies. They took 



a star formation rate per unit volume 



dp stars /dt = / H2 x (p sas /t SF ) x (ra H /[10 > 



,-3l\0.5 



(2.3) 



where igF = tg/eg is the free fall time, tg, divided by the fraction of the gas 
that turns into stars in a free fall time, eg — 0.02. To determine the molecular 
fraction, /h2, they considered heating and cooling, a radiation field proportional to 
the SFR, the ISternberg. et al. ( 2002) H2 formation theory, and radiative transfer 
using the Cloudy code. The result was a SFR that scaled steeply with the total 
gas column density, as observed, a higher KS slope for lower mass galaxies, which 
is also observed, and a shallower KS slope for the H2 column density alone, as 
in the molecular KS law. These results were somewhat independent of galaxy 



Bruce G. Elmegreen: Star Formation on Galactic Scales: Empirical Laws 11 



mass. The molecular /atomic ratio also scaled with pressure in an approximately 
linear fashion, regardless of galaxy mass, as observed. They also found a stability 
parameter Q that ranged from unstable in the inner, star-forming parts of the 
disk, to stable in the outer regions. 



2.2 Observations and Models of Outer Disks 



Murante et al. ( 201dh have a multi-phase SPH code that assumes pressure deter- 
mines the molecular abundance, and the molecules give the SFR. Below Etotai gas ~ 
10 Mq pc~ 2 , the slope of the molecular star formation law turns out to be very 
steep, Esfr oc £™ otal gas for n ~ 4. Above 10 M pc~ 2 , the slope i s the same as 
in the Kennicutt law, n — 1.4, which is steeper than in lBigiel et al. I (l2008h . where 
n ~ 1 for the molecu lar Schmidt law. 

Bush et al. ( 20101 ) simulated galactic star format ion with special atten t ion to 



the outer disks. The star formation model followed ISpringel fc Hernquistl (|2002h 
with radiative cooling, star formation in the cold phase, no specific molecular 
phase, and a volume- Schmidt law, psfr °c Pt t a i gas . They found patchy star 
formation in the outer parts, usually along spiral arcs where the g as density was 



high. This morphology is in agreement with GALEX observations (jThilker et al 



2005t Gil de Paz et al. . 12005 ). The Bigiel et al. and Leroy et al. observations were 



matched qualitatively in these outer parts too: below E g 



10 M Q pc" 2 , the 



slope n ~ 6 to 8 was steeper than in the observations (which also plot SFR versus 
total gas in the outer regions). Then it was les s steep at higher column density, 
with a slope of n ~ 1.4, which agrees with the Kennicutt] (|l998f ) slope for total 
gas. 

Outer disks can be Toomre-stable on average because the gas a nd star column 
densi ties are very low. This is especially true for dwarf galaxies ([Hunter et al. . 
1998). It might be that magnetic fields and viscosity destabilize outer disks, as 



discussed in Section l"!.!) but in any case, outer disks appear to be much more stable 
than inner disks. More importantly, the gas is outer disks is often far from uniform 
and the use of an average column density for Q is questionable. Locally there can 
be isla nds of high column d ensity where Q is small enough to be in the unstable 



region ([van Zee et al 



199G). These islands have to be larger than the Jeans mass, 
which might be 10 7 M@. Spir al arms and large disturbances in pressure could 
make unstable regions like this. Dong et al. ( 20081 ) found unstable islands of star 
formation in the outer part of M83. In these regions, the star formation followe d 



a steep KS law from point to point with a slope of about 1.4 ( Dong et al. . 2008) 



Bigiel et al. I (l2010al) found that outer disk star formation seen by GALEX fol- 



lows the HI very well in M83, with a uniform consumption time of 100 Gyr per atom 
beyond 1.5i?25- The form of the star formation is mostly in spiral arms. Outer 
disk a rms could be spiral waves radiating from the inner disk (jBertin fc Amorisco , 
l2010jL 



Boissier et al 



([20081 ) observed the galaxy-integrated KS law in low surface 
brightness galaxies, using a SFR from GALEX NUV observations. For a given total 
HI mass, the star formation rate was low by a factor of ~ 5 compared to normal 
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spirals, but over the whole range, the total star formation rate scaled directly with 
the HI mass. This is not the same as saying that the star formation rate per 
unit area scales directly with the HI column density because the observations are 
spread out in a plot like this with big galaxies on one side and small galaxies on 
the other. 



Bigiel et al.l (|2010bt ) studied SF rates in the far-outer disks of 17 spiral and 5 



dwarf galaxies, where the gas is highly HI dominated. The SF laws compare well 
with those in dwarf galaxies. There is no obvious Q threshold. They suggest that 
the total SF Law has three components, the extreme outer disk component that is 
HI dominated, a transition region where the molecular fraction increases to near 
unity, the molecular region inside of that, and the starburst component, where the 
surface density is higher than that of a single GMC. 



2.3 Scaling relations inside individual clouds 



Krumholz fc Tanl (|2007l ) showed that the conversion rate from gas to stars per unit 
free fall time is about constant inside clouds over a wide range of densities. This 
implies that the SFR per unit volume scales with the 1.5 power of density, with 
the first 1 in the power coming from the mass per unit volume, and the 0.5 in 
the power coming from the free fall rate. This is like a KS law, but for individual 
GMCs. There is also a thre shold column density for star formation inside GMCs 



of around ~ 5 — 7 mag in V (iJohnstone et al 
20061 IjQrgensen et al.l, koO 



2004: Kirk et al .. 2006; Enoc h et al 



Chen et al.l (|2010l ) studied the KS relation for individual GMCs in the LMC. 



They measured the star formation rate from both HII regions and by direct count- 
ing of young stellar objects. For YSO counting, the rate per unit area inside a 
cloud approximately satisfies the total-gas Kennicutt relation with the same time 
scale per atom, ~ 1 Gyr. For these regions, Shi+h 2 ~ 100 M pc -2 , larger than 
in the main parts of galaxy disks. Chen et al. also found that the areal rate of star 
formation was much lower in the long molecular ridge south of 30 Doradus than 
in the GMCs. Presumably this ridge is not strongly self gravitating, even though 
it is CO-rich. 



3 Summary 



The empirical star formation law on kpc scales is essentially one where star for- 
mation follows CO-emitting molecular gas with a constant rate per molecule, and 
the ratio of molecu l ar to atomic gas s cales nearly directly with the ISM pressure 
(jBigiel et al. , 2008; iLeroy et al.l . 120081 ) . The rate per molecule corresponds to a 
consumption time of molecular gas equal to about 2 Gyr. The place in a galaxy 
where the transition occurs between HI dominance in the outer part to H2 domi- 
nance in the inner part is at a pressure of P = 2.3 ± 1.5 x 10 4 fce K cm" 3 . There 
also tend to be characteristic gas and stellar column densities at this place, and 
a characteristic galactic orbit time for all of the galaxies observed. Beyond this 
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radius is the atomic-dominated outer disk. There, the SFR scales directly with 
Shi, and the consumption time is about 100 Gyr. 

Theoretical models of these empirical laws include the atomic-to-molecular 
transition in individual clouds and a sum over clouds to give the galactic scaling 
laws. Star formation occurs only in the densest parts of the clouds, as determined 
by a combination of turbulence-compression and self-gravity. Numerous simula- 
tions of star formation in galaxies can reproduce these empirical laws fairly well. 
The simulations usually show a sensitivity to the Toomre Q parameter, unlike the 
observations. 
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